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Abstract: We study the holographic dual model of quenched flavors immersed in a quark- 
gluon plasma with massless dynamical quarks in the Veneziano limit. This is modeled by 
embedding a probe D7 brane in a background where the backreaction of massless D7 branes 
has been taken into account. The background, and hence the effects, are perturbative in 
the Veneziano parameter Nf/Nc, therefore giving small shifts of all magnitudes like the con- 
stituent mass, the quark condensate, and several transport coefficients. We provide qualitative 
results for the effect of flavor degrees of freedom on the probes. For example, the meson melt- 
ing temperature is enhanced, while the screening length is diminished. The drag force is also 
enhanced. 
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1 Introduction and results 

Since the times of Rutherford's experiment, the most direct way of obtaining valuable infor- 
mation about the microscopic structure of a chunk of matter consists in sending hard probes 
that penetrate the system and, thereafter, analyzing the outcome. Nowadays, the quark-gluon 
plasma remains a challenging system for which a vast amount of data keeps accumulating 
along the last five years. The peculiarity of the heavy ion experiments performed at RHIC 
and CERN is that the probes are created spontaneously inside the plasma, and not introduced 
by hand by the experimenters. 

From these data, we have been able to see that many static and dynamical properties of 
probe quarks and mesons are dramatically aff'ected when immersed in a high density and tem- 
perature medium. Examples of phenomenologically relevant properties include the screening 
of the quark-antiquark potential, the melting of mesons (but not of glueballs) at a certain 
temperature, or the dramatic jet suppression eff'ects observed on back to back hard probes. 



- 1 - 



It makes sense to ask how much of these effects can be ascribed to the quarks that hve in the 
fundamental representation, and how much to the gluons, in the adjoint. 

The quark-gluon plasma should be described by a ground state of a theory which contains 
dynamical quarks and gluons. At least for temperatures at or above the deconfinement 
transition T ~ 170 MeV, the quarks u and d can be considered to be massless, and therefore 
fully dynamical. The strange quark, s, with mass ~ 10^ MeV should also be considered 
dynamical, albeit not massless in the same context. For the c and h quarks the assumption 
that they are quenched is justified by their heaviness, and the probes discussed in the present 
analysis may be interpreted as describing the phenomenology of these quarks. The t quark 
is too massive to be included in the regime of validity of our model, as discussed in the text 
(actually, the h quark is on the edge of the regime of validity). 

The problem of the propagation of probes inside a plasma is addressed in this paper by 
using the AdS/CFT correspondence. Our analysis is a continuation of the work initiated in 
[1]. The holographic setup consists of a solution to type IIB supergravity in the presence of 
D7 sources which are dual to massless matter transforming in the fundamental representation 
(quarks). It can be seen as a deformation of = 4 SYM where A'^^ D7 branes are introduced 
in the Veneziano limit: Nf — )• oo with Nf/Nc fixed. Using the presence of a large number of 
sources we smear them along their 2-dimensional transverse space, recovering some isometries 
of the internal space present in the unflavored type IIB solution. The geometry is dual to a 
field theory in which the beta function is positive. Therefore there is a Landau pole where 
a UV completion is needed. As we will review below, this is captured by the supergravity 
solution, and limits the current approach to phenomenology in the far IR limit of the theory. 
On top of this background we will consider a finite number N'j- <^ Nf,Nc of massive probe 
D7s. We are interested in the dynamics of these A^^ flavors, since they will give information 
about the behavior of massive quenched quarks in a plasma of dynamical massless quarks 
and gluons. 

The paper is structured as follows. In section 2 we shall review the background and its 
properties. We will also choose a criterium for comparing two theories, one without and 
one with dynamical flavor. The aim is to discriminate between replacing degrees of freedom 
(from adjoint to fundamental) from merely piling them up. In section 3 we will discuss the 
embedding of D7 probes in this background. The effect of the unquenched massless flavors is to 
induce logarithmic corrections to the profile which, in principle, should affect the expression 
for the {ipip) quark condensate. Deprived of a direct dictionary, in section 4 we will show 
how to extract the correct value of the condensate from the renormalized boundary action. 
Furthermore, by studying the condensate, we observed that the phase transition between a 
confined and a deconfined flavor (meson melting point) occurs at lower values of the quotient 
Mq/T. Finally, we also argue that the constituent mass can be computed reliably for masses 
not very large as compared to the temperature, and see how the presence of a sea of quarks 
increases it. 

In section 5 we address the problem of the conductivity of a U{1)b current on the D7 
probe. In order to do this, we shall accommodate for an additional non-vanishing value of 
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the ^0 component of the world-volume gauge field. This study suggests that the fundamental 
(quarkonic) degrees of freedom present in the background have larger cross section than the 
adjoint (gluonic) degrees of freedom. This point is supported also by the study of the drag 
force of heavy quarks. Furthermore, the pair production of quarks is shown to contribute less 
to the polarization of the vacuum than the pair production of gluons. 

Section 6 is devoted to the quark-antiquark potential. This is performed in the standard 
way by looking at the spatial Wilson line for long times, and confirms, as was found in similar 
cases, that the dynamical flavors tend to screen the color charge. 

Finally, in sections 7 and 8 we study the problem of energy loss of probes in the plasma, 
due to the drag force and the jet quenching effect respectively. In both cases the presence of 
fundamental matter enhances the loss of energy, thus suppressing back-to-back jets and stop- 
ping heavy quarks moving through the plasma. In this last case the heavy quark, modeled 
holographically by a string hanging from a brane near the boundary, feels an effective tem- 
perature which is larger in the presence of unquenched flavor with respect to the unflavored 
case, and this leads to a larger absorption of momentum by the presence of a larger black 
hole in the worldsheet of the string. 

2 The flavored A/" = 4 SYM gravitational dual 

The gravitational background we will consider was proposed in [1] as a holographic dual to a 
quark-gluon plasma with unquenched, massless flavors in the Veneziano limit Nc ^ oo and 
Nf/Nc finite, with Nc {Nf) the number of color (flavor) degrees of freedom. 

This scenario corresponds to a solution of type IIB supergravity in 10 dimensions with 
Nf D7 flavor branes. Those branes act as sources in the supergravity action for the Fi 
RR field strength, and backreact into the geometry created by a stack of coincident D3 
branes at the tip of a Calabi-Yau (CY) cone. The flavor branes are smeared along their 
transverse directions to recover (some of) the isometries of the Sasaki-Einstein (SE) manifold 
that defines the base of the CY cone. The construction is very general and can accommodate 
any 5 dimensional SE manifold X5, for which the metric can be written as a U{1) fiber bundle 
over a 4 dimensional Kahler-Einstein (KE) space 

dsx, = dsj^E + (dr + Ake? , (2.1) 

where Ake is the Kahler-Einstein potential, i.e., dA^E = 2 J with J the Kahler-Einstein form. 
Despite its generality, it is useful to keep two main examples in mind. The first one, = S^, 
is the usual original standard example of holographic dual to A/" = 4 supersymmetric gauge 
theory, where the base KE space is CP^. The second is the singular conifold = T^'^, 
which is dual to the Klebanov-Witten quiver, and the KE base is S'^ x S'^, in this case. 
The lOd metric ansatz introduced in [1] is given by the following line element 

ds^Q = Gttdt^ + G^^dxl + Grrdr'^ + G^^ds^^ + G^ridr + Ake? ■ (2.2) 
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The Fi RR field strength sourced by the smeared D7 brane squashes the fiber relative to 
the base in (2.1). Hence the difference G-rr — G^^ is a measure of the backreaction of the 
flavor branes onto the geometry, which will be parameterized by a dimensionless parameter 
eh- The configuration breaks the super symmetries from A/" = 4 to A/" = 1. Furthermore, finite 
temperature will break all supersymmetries. To keep the discussion short we will specialize 
to the case of AdS^ x 55. 

The solution to the sourced IIB system depends on a dimensionful (energy) scale^, r^, and 
a dimensionless parameter, e/,. The scale is the horizon location and sets the temperature 
scale T = rh/nR'^ + ©(e/^), with R the curvature of the manifold. The second parameter, e^, 
is proportional to the number density of fiavor branes 

where = AnggC^'^ N^. ^ 1 is the "running" 't Hooft parameter evaluated at the IR scale r/j. 

In the case of the five sphere, the SE manifold corresponds to CP^, so the KE metric is 
given by ds^g = ds^p2 + [dr + A)'^, where^ 

d4p2 = \dx^ + ^ cos^ ^{dO^ + sin^ Od^'^) + ^ sin^ xid^ + cosOd^f , (2.4) 
Ake = ]^cos'^^{di + cos ed^) . (2.5) 
For the ansatz (2.2), we have the perturbative solution [1] 



e/i , 2 5 1 rh 

k ft npr 



G,,=R'{l + r^ + ei^—--log^]+0{et)] , (2.8) 



G^ = «^^l-i+^i(i + ilog^J+0(4)j . (2.9) 

, = log (^) + (h (1 + 3 log ^} log ^ + ^ log. (1 - ^)) + O(4).(2.10) 

Here we absorbed the constant in the definition of the dilaton field (p = ^ — The 
Fi RR field strength has the form Fi = Qf{dT + Ake), with Qf a constant proportional to 
the number of "flavor" D7 branes, Nf. For the self-dual RR field strength one has the usual 
ansatz F^ = Qc(l+*)e5 with Qc a constant proportional to the number of "color" D3 branes, 
Nc, and £5 the volume form of the SE manifold X5. 



^Actually, the solution depends also on a UV scale, r*, acting as a Wilsonian renormalization scale, from 
which a completion of the theory is needed. We will comment later on this. 

^ The ranges of the angles in CP^ and the fiber are 0<x<i"i 0<^<7r, 0<(p<2TT,0<£,<'iiT 
and < r < 27r. 
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The above solution should be understood as an effective IR one in the Wilsonian sense. 
Namely, the perturbative solution that the reader can find in [1] contains another scale, r^, 
beyond which the solution itself should be replaced by the correct UV completion (r^, would 
be the analog of the Z boson mass Mz for the Fermi theory of electro- weak interactions) . One 
can show that the dilaton blows up at a radius rip ~ r*e^/'^''. From here we see that, besides 
the possibility of obtaining a perturbative solution, another reason to keep e/i ^ 1 is to have 
the scale of the Landau pole exponentially separated from the Wilsonian cutoff scale. The bad 
UV behavior of the solution is understandable on physical grounds, and signals the fact that 
the beta function is positive on the QFT side. Hence, the physical coupling constant is known 
to diverge at a certain UV (Landau) pole. In the spirit of effective theories, we expect that 
IR observables are mildly affected by whatever UV completion is added. In fact, as explained 
in [1], all corrections to physical results computed at a scale r < r^, are modulated by factors 
0{r/r^). Therefore, if we insist in investigating the far IR limit of the theory, r rh, we 
can discard the corrections coming from the UV completion. This has been realized formally, 
in the solution (2.6)-(2.10), by sending r^, — )• oo. The thermodynamic quantities that were 
evaluated in [1] are an example of IR quantities that do not depend on the UV completion. 
The temperature of this background is obtained from the absence of conical singularities after 
Wick-rotating the time direction by 

rh ( eh 13 2 



1_ - + . (2.11) 

Also the energy density can be adequately described within the present framework, giving 

'='-^^cT'{^+r^+ri+---) ■ (2-12) 

The present paper deals with comparisons of quantities as computed in two different theo- 
ries, one with flavor and the other one without flavor. Still there are several parameters that 
can be chosen to be equal in both theories. In [1, 2] an extensive discussion has been provided 
on the different possible choices. Here we will stick to one choice, by keeping the temperature 
T and the energy density e flxed. From (2.11) and (2.12) this simply amounts to shifting the 
value of r/i and Nc suitably. 

,(0) ^ ^9 .2 



JV, = N(«>(l-ii-^jJ + ...) , (Z14) 

with 

rr=vri?r, iVi°^ = ^/|, (2.15) 

the values in the flavorless background. In this way, rather than merely piling up degrees of 
freedom, we are partially replacing adjoints by fundamentals while keeping a measure of the 
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total sum of them unaltered^. For example, we will also need to shift the 't Hooft coupling 

A, = 47r5.e<^'^iV,(6,) = A^") (l - | - + • • • ) . (2.16) 

The fact that the solution (2.6)-(2.10) is perturbative in €h adds another source of divergence 
in the large r limit, coming from the truncation of the series at a finite order. Terms of the 
form e/j log(r/r/i) stop being perturbative corrections unless r <C ta = r^e^/'^^ (~ lO^r/j for 
Eh = 0.24 , see section 6 of [1]). In particular this implies we should not consider placing 
arbitrarily massive probe quarks in this background, as the point r^m of maximum approach 
of the D7 should be way below rmin ^ fh for the computations to be reliable. 



3 D7 probes on the flavored background 

In this section, and due to computational limitations, we will restrict ourselves to ther first 
order in for the backreacted solution. We will consider a class of embeddings for the probe 
D7 branes that, apart from extending along the Minkowski 4-space, wrap a non-compact 
cycle in the cone 

= To , X = X{r) ■ (3.1) 

For example, in the supersymmetric case, the embeddings are given by tq = rsin(x(r)/2), 
where the mass of the probe quark is proportional to tq and, hence, massless embeddings 
have X = 0- For a general profile, x(^)) the induced metric will read from (2.4) 

dsl = Gttdt^ + G^^dx^ + {Crr + ^x'^) dr^ + ^ cos^ | [d^^ + dO^ + dc/.^ + 2 cos 

^ Grr - ^Qg4 ^ cOsg#)^ . (3.2) 

Notice that, when the backreaction vanishes, Grr — Gj^^ = 0, the probe brane wraps an 
C , but this is no longer the case when e/i corrections are present. The lagrangian 
density for D7 probe branes in the Einstein frame is 

C = Cdbi + Cwz = Tr e* y/- det gs + Cs , (3.3) 

with Ty being the tension of the D7-brane and Cg is the RR eight-form potential (see below). 
Hatted are pull-back quantities in the worldvolume of the probe brane. Computing the 
determinant of the metric (3.2), the lagrangian density of the D7-brane reads 



Cos,= '^ sin cos' I ;-G„GLO,.OLl/ ' + (|j " f f + Ifx ' P'") 



^Alternatively, the entropy density could have been chosen but results remain unchanged up to first order 
in Efi, as s and e only differ at order e| 
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Let us now analyze the WZ term. First of all, we take the RR nine-form field strength Fg 
in the Einstein frame as Fg = — e^* -k Fi = —QfC^^ * (dr + A^e), where we substituted in Fi 
the form sourced by the background D7-branes. The corresponding Fg is 



Qf 



Fg = :iie2* sin 
^ 32 



' GttGrrGl^Gf^^ 



Gr 



cos^ ^ sin X d^x A dr A dx /\ d(f) A dO A . 



(3.5) 



We can define the potential Gg as -Fg = dGg and the corresponding WZ lagrangian density is 

(3.6) 



J-wz — H 



Qf 



,2* d3 , 



32 



GuGrrGl^G^ f 4 x\ 
^ (cos - I sm ( 



Rewriting the dilaton by expressing explicitly its value at the horizon, e*^^) = e*''"'"'^^''^ and 
using the definition e/j = Q/e^'^, we arrive at the total 8 dimensional lagrangian density 



C 



-e^ sm w ^|Gt^GrrG\^G\ 



(3.7) 



X COS — 

2 



'1 + 



X 



cos — H — - e 



G-T-j- 2 



It will be useful to make the following change of variable ip{r) = sin(x(r)/2) in order to 
make contact with the notation used in the literature^. The dimensionless action density, 
Id7 = S/VolCRs), acquires the following form after integrating in t G (0,/3 = 1/T), and 
6, (j) and ^ 



lD7=M{eh) / dr^JGuGrrGl^GlJ 



1 + 

Gr"r 



1 + 



1 



4 



G. 



(3.8) 



where we have defined the short hand notation -0 = ^/l — ip'^. Notice than in (3.8) the integral 
is dimensionless thanks to the factor r^^ which has been added by hand for convenience. The 
prefactor, which gives the right dimension 3 to the action density Idt^ reads 



with 



^<«.(i,f^,JL,j,...). 



(3.9) 



(3.10) 



*This is exactly the same function that is called x = cos 5 that appears in [3]. Clearly massless embeddings 
have %p = 0. 
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Expressing (3.9) in terms of field theoretical quantities was made by means of the relations 
Tfi = {27TygJl i?4 = ATTgsNJi = g^^NJ^ and A;, = Airg^e'^'^N,, as well as (2.13)-(2.16) 
which will be used extensively during this note. 

The equations of motion for the embedding can be derived from (3.8). After inserting the 
perturbative solutions given above for the metric components and dilaton, the equations of 
motion can be expanded in powers of e/j. A perturbative solution for the embedding involves 
an ansatz 

oo 

Hr)=Y.<Mr) , (3.11) 
p=0 

where each function ipp{r) solves the equations of motion at a given order e^. As mentioned 
above, we shall only compute up to V'i(^) ^-iid consider this as the first order correction to 
the physics of the probe induced by the fiavored background. 

4 Flavor corrections to mass and condensate 

In the previous section we have displayed an effective gravitational solution where we have 
decoupled the UV Landau pole by sending it to r^. = oo. In [1] this solution was shown to 
provide sensible results for IR quantities of the dual field theory. In particular the thermo- 
dynamics as well as the hydrodynamics [4] are found. It seems natural to expect difficulties 
when addressing questions which involve a proper treatment of the UV divergences through 
holographic renormalization. It is the purpose of this section to elaborate on such aspects for 
the probe branes. 

Central to the discussion is the choice of integration constants. The perturbative expansion 
of the equation of motion for the embedding profile, il^{r), leads to a tower of second-order 
ordinary differential equations for the components 'ijjp{r), p = 0, 1, • • • . A series solution 
around the UV, r ~ oo, depends on two constants, related to the bare mass of the quark Mq, 
and the condensate {ip^) in a way that we will clarify. The value of the constituent quark mass 
is a natural boundary condition to choose and, naively, the value of the condensate would 
be the second condition to impose. However, a generic value of {ipip) leads to a divergent 
solution for the bulk field ipp{r) in the IR of the theory. This can be avoided by imposing 
regularity in the IR as the second boundary condition. 

Numerically it is easier to start the integration from the IR towards the UV. Then regularity 
is satisfied by construction, and the series solution around the IR limit depends upon a single 
integration constant, ^, which parametrizes the condition we impose there. For example, 
consider the unflavored component ipo{r). Depending on whether we are dealing with a black 
hole or a Minkowski embeddings [5], the integration constant can be the insertion angle of the 
D7 brane when it pierces the horizon ^ = ipo{rh), or the minimum distance ^ = rmin between 
the probe D7 and the color D3 branes, defined by ^o(^mm) = 1, respectively. 

Concerning ipi{r), the first order correction in to the embedding profile, it is governed 
by a second-order, linear, non-homogeneous differential equation. As before, we integrate 
numerically from the IR towards the UV with adequate boundary conditions. For black hole 
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embeddings we chose ipi{rh) = and for Minkowski embeddings V'i(^mm) = 0, so that both 
the anchoring point and the minimal distance are fuhy encoded in V'o(^)- 

Up to the first correction in flavor backreaction the embedding profile has an expansion for 
large r given by 

i/jir) = i^oir) + ehMr) + ■ ■ ■ (4.1) 

rh \ , frh\'^ , / mo . r frh\'^ f 5co , r\\ 

= — mo + — Co + eh [ mi - — log h — ci - — log — + • • 

r I \ r / \ rh ^ r / \ b rh J J 

where the constants nio, mi, cq and ci all depend on ^ and can be extracted from a fit to 
the numerical integration. We see the appearance of logarithmic terms in the expansion. A 
similar effect has been observed in various situations before. In [6, 7] it was seen to arise from 
a coupling to the Ricci scalar in four dimensional space. Our situation here is more akin to 
the one in [8]. Observe that we obtain logarithms contributing both to the mass term and to 
the condensate, so the proper identification of lagrangian mass and condensate are suspect. 
However, we shall argue later that, for such an embedding, the quark bare mass is indeed 
given by the sum^ m = mo + e/i"^i- 

Conversely, we may be interested in seeing how the D7 profile is modified by the addition of 
flavors for a fixed value of m. This amounts to shifting £, ^ £,' = £, — e/i^i, in such a way that 
m{^') = mo(^) stays invariant. This is important, for example, in order to deduce how the 
presence of flavors modifies the constituent mass, Mc ^. This is only defined for Minkowski 
embeddings, and is obtained from the Nambu-Goto action for a fundamental string that 
stretches from the horizon up to the point of nearest approximation of the D7 brane to the 
origin. The relevant calculation is 



1 r'^min 



eW2 
2W 



-, 3£/i \ , \ , ^h 1 r^Yiiji 

1 — ] [rmin -rh) + —rmin log \ 

8 / 2 Th 



(4.2) 



In the flavorless limit eh = 0, the above expression allows for a physical definition of the bare 
quark mass. In that case, the UV regime is safe and conformal and, thus, in the limit of very 
high masses as compared with the plasma temperature T, the constituent mass Mc and the 
quark mass Mq should become the same. From (4.1) in the unflavored limit, we see that the 
definition of the minimum distance leads to r^m = mrh + 0{l/r). Hence 

This result also holds in the limit T — >• 0, as can be seen from the exact super symmetric 
embedding: mrf^ = rip{r). 



^In fact, one can perform a reparameterization of the radial coordinate r ^ f = r{l + 0{eh) + ■ ■ ■ )i such 
that the logarithmic term that shifts mi disappears. 

®In [9] the constituent mass is denoted rest mass Mrest- 
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When Eh is not zero this nice matching is lost. Naively we would expect the large mass 

limit to give Mc — ?• \\J xf^^T{mQ + ehTni), since a very heavy quark should decouple both from 
thermal fluctuations as well as from the massless quarks in the sea. However taking the limit 
of large rmin in (4.2) gives a divergent result. This is, as adverted, one of the drawbacks that 
arise from truncating the perturbative expansion in powers of e^. Hence M(.{eh) cannot be 
trusted beyond rmin ~ = r^e^/''^. For an eh ~ 0.24 (which would correspond to a realistic 
regime with ~ 67r and Nf = Nc = 3) this gives a safe range for rmin/i'h ~ ^ 0(10). In 
this range we must anyway resort to numerical analysis. 

M 




Figure 1. In this figure we plot M/{^y/\T) for different M's. For the black diagonal line M = Mq 
the lagrangian quark mass. In blue we see M = the constituent mass in the unflavored thermal 

background. It vanishes for rmin ~ 0.92. The purple curve is for M = Mc, the constituent mass in the 
flavored thermal background. It seems to cross the black line at a finite value of m. This is an artifact 
of the approximation and does not happen for the restricted range of m. In a realistic model with a 
controllable UV completion both blue and purple curves should approach the black line for m — >■ oo 



This means in the context of the present analysis, we should only consider quark masses 
bound as Mq = \^/\^Tm < 0(20 T) ~ 3.5 GeV (taking for T the critical temperature 
Tc = 175 MeV as the lowest possible physical value). In the plot of figure 1 we see the 
influence of the massless unquenched flavors on Mc of the quenched probe in the advocated 
restricted range. The generic lesson is that fundamental degrees of freedom contribute more 
than adjoint ones to the constituent mass. 

4.1 Holographic renormalization, quark mass and condensate 

In this section we will justify the aforementioned identiflcation of the bare mass. We will 
moreover provide an expression for the quark condensate. First of all, let us deal with the 
on-shell action. After inserting the asymptotic expansion (4.1) into the action (3.8), we 
obtain quartic and quadratic divergences that can be tamed with a standard counterterm 
flrst proposed in [10] and, in a similar context to ours, in [11] 



/c,i=AA(e,)v^e^|^(l + |) {l-Hr?y 



(4.4) 
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with 7 the determinant of the induced metric on the r = rc hypersurface, which is a convenient 
(large) cutoff. The subtraction Ijjj — Ict^i removes ah power hke divergences, but leaves a 
logarithmic one. This last can be eliminated with an additional counterterm of the form 



ct,2 



^^/i log 



rh 



(4.5) 



It is remarkable that this divergence is proportional to rh and, hence, only pops up upon 
turning on the temperature. In the supersymmetric case, the addition of (4.4) is just enough. 
It is unclear to us what can be the source of this new divergence, although there might be 
several possibilities (see for example [12, 13]). The important fact to stress is that, unlike the 
case in [8], this counterterm is totally independent of the details of the probe brane profile. 
Therefore, it will not affect the computation of the condensate. 

Let us now proceed to compute the 1-point function. One would naively consider (ipip) = 
Tdlren/dm with m = mo + emi as a correct prescription for the condensate. In this sense 
the renormalized action Iren = hulk ~ Ict,i is not useful, as it does not exhibit the form of a 
product of source times v.e.v. ~ J J{0) where J ~ mo + euf^i- As very well explained in [3], 
the point is that both cq and ci depend implicitly on mo and mi through the fact that all 
of them are controlled by a single parameter: either ^ = rmin for Minkowski, ox ^ = iph for 
black-hole embeddings. Hence it will be convenient to think of Iren as parameterized by ^. 
We compute I'reniOi yielding 



I'renii) = I'niiO - I'etM) 



dC di) 



diet A dip 



dip 



(4.6) 



where rc is a cutoff radius, use of the Euler Lagrange equations of motion has been made, 
and an integration over all coordinates but r has been understood in C. As said before, Ict,2 
does not contribute to this calculation since it does have any dependence on the embedding 
profile. The contribution at rmin vanishes, and after adding up both terms the final result 
organizes in a neat form 



2cq + eh ( 2ci + ^cq 



{m', + ehm[) + 



(4.7) 



Using the chain rule this gives the holographic dictionary we are after, if we interpret the 
bare quark mass as Alq = ^ v/ X^j^^Tm with m = (mo + €hmi) 



u±ren 



/ 



2co + eh ( 2ci + -Co 



(4.; 



In figure 2 we have plotted the values of the condensate versus the mass. It is intriguing to 
observe that the effect of the massless backreacted flavors is to lower the curve in such a way 
that the condensate changes sign at a given value of the mass, whose position depends on eh- 
A similar effect occurs upon turning on a magnetic field B as found in [14, 15]. Unlike those 
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Figure 2. Here we observe that the effect of the flavor mass density on the heavy quark condensate 
is to lower its value. From top to bottom the curves correspond to values of = 0,0.2 and 0.4 
respectively. On the right hand .side we take a closer look at the transition region and observe that the 
critical mass wicr is reduced, as can be deduced from the Maxwell equal-area construction 



papers, here we cannot see if a spontaneous chiral symmetry breaking occurs for high values 
of efi, as this parameter is perturbative in our approach. The presence of unquenched matter 
also shifts the point of the meson-melting transition, such that for fixed Mg the transition 
happens at a higher temperature^. 

5 Conductivity and transport 

The presence of massless fundamental degrees of freedom in the background has an impact 
on the transport properties of the probe brane. Fluctuations of the fields defined in the 
worldvolume of the probe D7 are sensitive to the Ch corrections in the background fields 
(i.e., metric and dilaton), and so is the response of the system to these perturbations. The 
relations between these responses and the fluctuations sourcing them define the transport 
coefficients. A paradigmatic example is the (microscopic) conductivity, which relates, with 
the aid of linear response theory, the microcurrent induced in a system to the fluctuation of 
the electric field that causes it. In a homogeneous and isotropic medium this relation reads 
— (^microi^)^E [uj) , where uj is the frequency of the fluctuation. 
There are examples in the literature where the calculation of (Tmicro has been performed 
in different backgrounds. We refer to (some of) those works for details of the calculation 
that we only sketch in the this note. We will restrict our study to the DC limit of the 
conductivity an alternative approach can be used, as devised in [16], with a clear physical 
significance. This approach implies a macroscopic calculation of the conductivity in which 
one mimics an Ohm's law experiment. Instead of studying fluctuations of the electric field in a 

''Notice that, as compared with [3], our A — ^tvQs is a factor of 2 larger. Expression (4.3) however is the 
same because in (4.2) we have absorbed a factor of 1/V2 in mo, mi, as well as a factor 1/(272) in Co , ci . This 
explains the discrepancy in the meson melting point, which is at mo — 0.923 for us instead of being at 1.305 
for them. 
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setup where no background E is present, one considers a worldvolume gauge field of the form 
A = AQ{r)dt + {Ext + Ax{r)) dx . In this expression E^ is a constant electric field, that we 
align with the x direction. The radial dependent fields Aq and Ax determine holographically 
the baryon density, -/ig, and charge current, (Jx), in the field theory [16]. From Ohm's law we 
define the macroscopic conductivity as {Jx) = (TmacroEx ■ Using the membrane paradigm it 
has been shown that the macroscopic and microscopic conductivities are related as 

hm (Jmacro{Ex) = hm Gmicroioj) = (TdC , (5-1) 

but the macroscopic calculation is much neater and direct. 

The calculation can be done for black hole embeddings only, since Minkowski embeddings 
represent an insulating phase where the conductivity vanishes due to the discrete spectrum. 
Following the standard procedure described in [16], we obtained the following result: 



<^Dc = agi - 1) (1 - i^i? + (i - 1 j ' (5-2) 



where V'/i = '^{I'h) and we have defined 



-St = , = . (5.3) 



To plot expression (5.2) we do not need to calculate the embedding profile in the presence 
of finite baryon density, if we want to plot the conductivity as a function of iph- This step 
becomes mandatory, though, if we want to trade the IR quantity, -i/^/i, for the physical UV 
one, m, as the independent variable in fig. 3.® 

As pointed out in [16], two different mechanisms for the conductivity can be identified in 
(5.2). One of them comes from the presence of a density of free charge carriers proportional to 
the density charge. Its contribution gives rise to the d^^^ term under the square root. The other 
one remains even in the charge-less limit d^^^ — t- and is proportional to the temperature- 
dependent factor o"g^, signaling a thermal origin. The interesting fact pointed out in [16] is 



* Obtaining the modified embedding profile goes like in section 3, but with a non trivial gauge field 
A = Ao{r)dt. The presence of this field in the D7 action shows up in two places. The first one is in the DBI 
part of the action, which expands from ^— det g into \J~ det(,g + 2-Ka' F). As a consequence, the last square 
root in equation (3.4) should be modified as follows 



1 I X" . L , G^, x'^ {2na'YA',^ ^^^^^ 

Grr 4 y Grr 4 \Gtt\Grr 

The second contribution has to be considered when the unquenched fundamental matter has a net charge [2]. 
Then, in the Wess-Zumino part of the probe action there is a term coming from C& A {A'^dt A dr), where Ce 
is the puUback of the background Ce RR potential to the worldvolume of the probe brane. This gives a term 
proportional to JA'q, with J the function taking into account the effects of the backreaction of charged D7 
branes (see reference [2] for more details on this function). We will stick to the case in which the background 
is not charged, therefore J = and the WZ contribution does not get modifications. 
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Figure 3. ctdc/'^^dc '^•^ function of the (renormalized) constituent quark mass parameter Mq for 
values of e ~ 0, 0.2, 0.4 (from blue to red) and d ~ 0.01, 0.7 and 1.2 from bottom to top. 



the presence of the modulating factor (1 — V'h)^^^ which suppresses this contribution for high 
masses {iph !)• This points to a Schwinger-like pair production of charged carriers as the 
mechanism behind this effect. 

In (5.2) the correction to these terms due to the backreaction of flavor degrees of freedom 
is a decreasing one. For the term depending explicitly in d^^\ the electric field accelerates the 
valence charge carriers, and the drift of charged matter gives rise to a macroscopic current. 
As they move, they scatter against the bulk adjoint and fundamental degrees of freedom that 
can absorb momentum and energy from these carriers and, therefore, slow them down. The 
negative correction (1 — eft/4) means that the fundamental degrees of freedom have a larger 
cross section than the gluonic ones they have replaced. As we will see, this is in agreement 
with the enhancement of the drag force on a heavy test quark. 

Concerning the first term, we see that the correction is also negative (1 — e^/G), whereas 
we would expect the pair creation to be enhanced by the presence of charged massless fun- 
damental matter that can run inside the loops. Here we have to stress the importance of 
the comparison scheme selected. We have chosen to keep invariant the energy density, and 
therefore the number of colors is diminished according to (2.14). Had we decided instead to 
compare two theories with the same number of colors, and just additional fundamentals, the 
correction to the first term would have been (1 + e^/S), hence positive in agreement with the 
aforementioned expectations. From the fact that the net correction is negative, we conclude 
that at fixed number of degrees of freedom the fundamental matter contribute less to the 
pair-production than the adjoint matter. 

In figure 3 we plot the conductivity as a function of the quark mass m = {niQ + ehmi) for 



- 14 - 



various values of the density d^^\ For d^^^ — )■ the curves drop to zero conduction (insulator) 
at a point that drifts from m ~ 0.92 downwards with Ch- This is consistent with the shift in 
the meson melting transition point observed in fig. 2. We observe that at large mass the pair 
creation term is negligible and the conductivity approaches a constant d{l — e/8). 



6 Quark potential 



It is well known that one can extract the potential energy of an external quark-antiquark pair 
from the expectation value of a Wilson loop operator 



{W{C)) = yl(L)e-^^(^) 



(6.1) 



for a rectangular loop which is an infinite strip T — t- oo. Here L is the quark-antiquark 
separation in Minkowski space. The classical calculation in [17-20] can be easily adapted to 
the present setup. Since the only dimensionful parameter is the temperature, we expect on 
general grounds a screened Coulomb potential of the following form 



V{L) = (1 + c(Lr)4 + e/(LT) + •••), 



(6.2) 



where c is a constant [20]. The idea here is to parametrize both L and V through the hanging 
string maximum depth ro 



L 



dr 



>"min(<:h,m) 



dr 



R\/r^-r1: 



1 + ^(1-4 



4 4 



I log- 



ro 



1 + 1 + 4 loo 



r 



4 log 



+ 



r 



, (6.3) 



(6.4) 



1 + ^ ( l + 41og 



r 



+ ■■■ 



In this expression, rmini^h^rn) is given by the point of minimum radius for the probe D7 for 
a given value of the bare mass parameter m. It has to be obtained numerically by shooting 
from different minimal values of rmin until, in the UV, the desired value of m = mo + e^mi 
is hit. 

In figure 4 we show the curves that correspond to the gg-potential for three values of 
e/j. The curves have typical shape obtained in [19, 20], which have been interpreted as the 
thermal screening of the plasma over the Coulomb potential. They exhibit a cusp at a certain 
separation L^ax, beyond which the U-shaped string embedding ceases to exist. This length 
is however larger than the dissociation or screening length Lgc, were the curve E{L) changes 
sign. Beyond this point, the dominant saddle point is that of a pair of disconnected strings 
hanging vertically down to the horizon. 
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Figure 4. The quark- antiquark potential Eqq as a function of the interquark separation L. The 
usual plot shows a turning cusp point where is ceases to exist. Before that, the configuration with 
two disconnected strings dominates whenever the energy changes sign. This defines a screening length 
Lsc = 0.37, beyond which the interquark potential flattens. As we see in the figure, the main effect of 
the flavors is to decrease this length. The plots are for eh = 0, 0.25 and 0.5 from blue to red. 

In the figure, we observe how the addition of massless flavors moves the curve towards the 
left. Hence the onset of a flat potential occurs at a closer separation that in the unflavored 
case. This behavior is in agreement with the result obtained in [21] for the same calculation 
in the Klebanov-Witten model at zero temperature. We also notice the truncation of the 
curves in the left plot of their fig. 3, due to the horizon which triggers the screening by 
string splitting (the saddle point with two hanging strings becomes favorable). Hence the 
usual double valuedness leading to the cusp is in the unphysical region. This is a direct 
consequence of the fact that, keeping the temperature fixed, demands rising the value of the 
horizon radius, r^, in (2.13). 

7 Drag force 

Let us now turn our attention to the issue of energy loss of partons in this medium. In the real 
world, the QCD quark-gluon plasma appears to be strongly coupled [22]. Moreover a very 
efficient mechanism of energy loss must be at work in order to explain the phenomenon of jet 
quenching [23]. There are two main ways of addressing the computation of this coefficient 
within the holographic paradigm. One makes use of an eikonal approximation at high energy 
yielding a non-perturbative definition of q as the coefficient of in an almost light-like 
Wilson loop with dimensions L~ ^ L. The implementation of this proposal in the string 
theoretical framework was done for the first time in [24]. 
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The other involves the picture of a quark travehng through the medium modeled entirely 
within a string theoretic framework. A parton moving inside the plasma is described by a 
macroscopic string attached to a probe flavor brane placed in the background that contains 
a black hole. The interaction with the medium is encoded in the string profile that dangles 
down the bulk and extends infinitely along the horizon. One searches for solutions for which 
the dynamics of the string endpoint, hence the quark, can be compared with a solution to 
the drag force equation 

P = -fJ'P + f, (7.1) 

where fi is the friction coefhcient, related to the mean decay time /i = 1/r, and / is an 
external force on the quark. As very neatly explained in [9], there are two simple solutions 
that one can readily device. The first one is stationary p = as a consequence of the constant 
worldsheet electric field that pulls the string endpoint 



X^=vt + xir), X' 



0, 



2,3 . 



(7.2) 



The string description of the dragging phenomenon embodies interesting predictions. Among 
them, an intriguing one is the appearance of a black hole in the pullback metric, for finite 
velocities of the quark. The horizon sits at a different radial position than that of the bulk 
metric, = As a consequence of this, the worldsheet and bulk Hawking temperatures 

differ by a velocity dependent factor Ts{v) = T/^ + 0[eh). This worldsheet black hole is 
at the origin of the thermal fiuctuations that are responsible, in the present setup, for the 
transverse momentum broadening of the heavy probes. Moreover, as remarked in [25, 26], this 
is the temperature that should be used in the relativistic version of the Einstein relations^. 
A general expression for arbitrary backgrounds can be found in equation (B.9). The flavor 
modification factor for this quantity is as follows: 



V7V 8 384 ^ 



(7.3) 



The worldsheet temperature in the flavored background is higher than the supersymmetric 
limit (but still lower than the bulk temperature, as it is expected and shown in flgure 5). 

To keep the velocity constant, the momentum absorbed by the charged quark, sitting at 
the string end-point, is dissipated into the bulk and through the horizon by the trailing part 
of the string, at a constant rate. It is easy to work out a general expression for the drag force 
in an arbitrary (dilatonic) background (see also [27] and [26]). 



/ 



2'Ka' 



\l Gxx I Gtt I 



(7.4) 



^As we will explain in the next section, we checked that Ts is indeed the temperature appearing in the 
Einstein relation for the correlators, even in the backreacted background, by using the Schwinger-Keldish 
formalism. 
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Figure 5. On the left we plot the worldsheet temperature as the quark moves inside the flavored 
plasma, relative to the unflavored case. On the right we see the modification of this worldsheet temper- 
ature over the bulk temperature is suppressed in the region v ^ 1 by the 7 factor. The different curves 
correspond to increasing values of = 0,0.25,0.5,0.75, from blue towards red. 





Figure 6. On the left (right), the ratio of the drag force (diffusion constant) with respect to the value 
it takes in the unflavored background is plotted as a function of the velocity, for different values of the 
perturbative parameter th = 0, 0.25, 0.5, 0.75. 



As will be the case for quantities that do not need the D7 brane profile, we can give a result 
up to second order in e/, 

^2 



/ = (1 + |(1 - log(l - v^)) + ^ (15 - 141og(l - v^) 

+91og2(l-t;2) + 121og2f^) ) +••• , 



(7.5) 



with 7 = l/Vl — v"^. Postulating a relativistic dispersion relation of the form p = Mj^in^v 
for some kinetic mass Mi^i^i and inserting in (7.1) (with p = for a stationary solution), the 
above calculation of the drag force implies knowledge of the joint product of fx and M^j^ 



/ 

7t> 



(7.6) 
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In the unflavored case, dual to A/" = 4 supersymmetric Yang Mills theory, a value of fj,^^^ was 
obtained that was independent both of the mass and of the velocity of the quark 

^(°)MS = |V^T^ (7.7) 

a fact that contradicts the weak coupling prediction [28]. Prom (7.5) and (7.7) we see that, in 
the flavored background, fj, acquires a logarithmic dependence on the velocity. The product 
fiMkin is related to the quark diffusion constant: 

D = , (7.8) 

where, as mentioned above, Tg is the correct temperature to be used here. Of course this 
is only true if the Langevin equation is a realistic description of the difusion phenomenon. 
In fig. 6 we see the modification of this constant with the flavor density e/,,. In summary, 
fundamental flavors contribute more than adjoint matter to the drag force and they inhibit 
more the diffusion 



D 



3 - 47^ + 107^ + (67^ - 87^) log(l - v^) + 37^ log^(l - v"^) + 127^ log2 v"^ ] + ■ ■ ■ 



The friction coefficient fi is related to the relaxation time r = 1//U. In order to disentangle 
fi from Mkin the authors of [9] propose to consider a different kind of experiment, in which 
the string momentum decays exponentially to zero in the absence of drag. Then / = in 
(7.1) implies that p{t) = p(0)e~^* for the probe quark. In the string picture this involves a 
quasinormal mode computation: the string has to fulfill outgoing boundary conditions at the 
horizon, and Neumann ones at the D7 brane. This cannot be done analytically for generic 
masses and, by construction, it only yields // in the slow v ^ limit. 

By solving analytically the limit of small masses Mc — )• (m — )• 0.92), the authors of [9] 
postulated the existence of an upper bound for the friction coefficient 

< 2ttT . (7.10) 

For nonzero flavor density, we observe from flgure 7, that the friction coefficient is decreased 
and, in particular, the bound is respected. 

Combining the knowledge of // with that of /iM^j^ previously calculated in (7.6) allows 
to solve for the kinetic mass. This is, by definition, velocity independent, and should only 
depend upon the bare mass m. The result is shown in the plot 8. Both the constituent 
mass, and the kinetic mass, increase as the fiavor density is switched on. The blue curves 
were obtained in [9] and the continuous (dashed) blue curve represents the kinetic mass M^j^ 
(constituent mass Mi°)) at zero flavor density. M^^^ is a monotonous function of m which 
has a minimum, for m 0.92 at VAWr/4. The purple curves plot the same magnitudes for 
e/j = 0.25. These curves are only trustable in this interval of masses. 
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Figure 7. On the left we show the dimensionless quantity /i/TrT as a function of the quark bare mass. 
The plots stop at m ^ 0.92 where the constituent mass vanishes. The horizontal line indicates the 
conjectured upper bound for this quantity. The ratio of the drag coefficient in the flavored background 
with respect to the conformal value is plotted on the right, as a function of the quark bare mass m. 
The different colors correspond to different values of the parameter Ch = 0, 0.25, 0.5 (from blue to red). 




Figure 8. The plot on the left shows the constituent mass Mc (dashed line) and the kinetic mass 
Mkin (plain line) as a function of the bare mass m. On the right we represent the ratio of the kinetic 
mass in the flavored background with respect to the kinetic mass in the unflavored background, as a 
function of m. Different colors correspond to values of the parameter = (blue), 0.25 (purple) and 
0.5 (red). 

8 Jet quenching 

The jet quenching parameter measures the momentum broadening of a quark traversing the 
quark-gluon plasma. If we assume that the heavy quark undergoes a Langevin diffusion 
process inside the plasma, the diffusion constant associated to this process is simply related 
to the jet quenching parameter. At a classical level, the Langevin diffusion reduces to an 
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ordinary diffusion process. 

The Langevin quantum description implies the presence of a stochastic force acting on the 
quark, both in the direction longitudinal and transverse to the motion of the quark. Let us 
denote the two components of this force by respectively. The two-point correlator of 

this force defines the diffusion constant k. Namely, at sufficiently long times (with respect to 
the typical correlation time of the medium, roughly the inverse temperature), we can write 
{F{0)F{t)) = K6{t), for each longitudinal and transverse direction. 

As it has been shown [29], this allows to determine the momentum broadening of the quark 
in terms of the diffusion constant. We adopt the definition of the jet quenching parameter as 
the average momentum squared of the quark divided by the distance the quark has travelled, 
q = ijP') /L. Keeping in mind the definition of the diffusion constant, the average momentum 
squared in the longitudinal and transverse directions can be obtained in terms of Ky and 
K_\_ respectively. This means that, for large enough times t, the quark acquires a transverse 
momentum that averages to zero but has quantum fluctuations (p^) = 2K±t. The longitudinal 
momentum has a classical part related to the drag force and quantum fluctuations from the 
stochastic force (Ap|) = Kyt. 

The diffusion constants can be computed holographically using the prescription for the 
holographic retarded correlators [30]. In fact, they are given by the zero frequency limit of 
the imaginary part of the retarded correlator, k = 2Tslim^^oI™ Gr{<^)- The procedure 
implies finding the solution for the trailing string fiuctuations with appropriate boundary 
conditions and evaluate the bulk-to-boundary correlator that, in turn, yields the retarded 
correlator. We apply the steps of [25, 29, 31], where the authors obtained the diffusion 
constants for the = 4 plasma, to the background with backreacted flavors. The aim of 
this section is to compute the jet quenching parameters in the flavored background using the 
following formulae: 

9± = 2— , q\\ = . 8.1 

V " V 

In order to compute the retarded correlator corresponding from the trailing string, the 
string ansatz must include the quantum fluctuations 6X^ 

= vt + x{r) + 6X^{r,t), X" = (5X"(r, t), a = 2, 3. (8.2) 

By a change in the time coordinate t = t+(^{r), with C'('') = dtr/stt [32], we obtain a quadratic 
action for the fluctuations with no mixed dtSX drSX terms 

•^^'^ = "2^ /rft dr ^ [gr{dr6X^f + gr{dt6X^f] , (8.3) 

1=1 

with 



(-G^,,H-^ \ 
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and 



Tj ± I Gxx {GttGxx — '^Gtt,sGxx,s) {Gtt,sGxx — GttGxx,s) /„ r-\ 

n-e^i „ „ „ , 



^ ^ ^iii^Gxx^ GttGxx - e'f"'-'t'Gu,sGxx,s ^^^^ 
Gxx y Gtt,sGxx ~ GttGxx, s 

Here the subscript s indicates functions evaluated at the worldsheet horizon rg. The retarded 
correlator is then obtained as 

Gr{u;) = - [^&Jj(r,a;)g-a,,M/«(r,u;)]b,,„dary > (8-7) 

where \I'/{(r, cj) is the retarded solution to the fluctuation equation, ^(r, oj) = e'^'^^5X{r,i), 
with boundary conditions given by ^/j(r, a;) — )• 1, at the boundary, and ^ij(r, w) (r — 
^^yui/i-n-Ts ^ at the worldsheet horizon. The fact that the fluctuations behave like infalling waves 
with temperature Tg at the worldsheet horizon implies that Ts appears in the aforementioned 
Einstein relation. 

In fixing the boundary normalization condition, we assume that the fluctuation has a 
regular expansion in the UV. By this we mean that for large radii the leading behavior of 
^ is constant. Indeed, using the expressions (2.9)-(2.10) and expanding for a large enough 
r such that ^ r <C r/jC^/'^'', we arrive to the expected expansion for the wave function: 
^ — C'source + Cvcv&~'^^'^ / ■ The normalization condition can be set at a cutoff radius r^m 
(related to the probe quark mass Mc by equation (4.2)). This avoids the problem of going 
beyond the truncated expansion in Ch for the bulk fields and we can stay in the regime where 
T^s ^ T^min ^ rhB^^^'^ . Moreover, the zero frequency limit of Im Gr is independent of this 
cutoff. Hence we can consistently derive the diffusion constants (and consequently the jet 
quenching parameters) with a cutoff and arrive to a final cutoff-independent result. 

In order to compute the zero frequency limit of the imaginary part of the retarded cor- 
relator we do not need the full solution, but only its asymptotics. Moreover, since Im Gr 
is a conserved flux, it can be evaluated at an arbitrary point in the radial coordinate (not 
necessarily at the boundary). In particular, it is easier to evaluate it at the horizon, using the 
in-falling wave function expression ^R{r) ~ ^'^(r — r^)*'^/^'^"^" . Substituting this asymptotic 
into the formula (8.7) for the retarded correlator and reminding the definition of the diffusion 
constant, we obtain 

1 



jBX\^s\^ (8.8) 
Tra 

with Bs defined by 

g""(r)~47rT,S2(r-r,) + 0(r-r,)2, r ^ r,. (8.9) 

The constant is evaluated by matching the low frequency limit of the retarded solution with 
the exact zero frequency solution and imposing the normalization condition. Summarizing the 
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results of this computation, we get ^'^ = 1. Moreover, deriving the asymptotics of equation 
(8.4), we obtain i?^ = e'^/'^G^x and similarly for = e'^'/'^G^xZ'^ ■ Therefore, the diffusion 
constants and jet quenching parameters read 



<1± 



2^T^^/W)^ (l + ^(l + ^2^21og7)+ 



^.10) 



1\\ 



.2 \ 

+-^(15 + 14^2 - + 4(7 + 3t;2) log 7 + 36(log 7)^ + 12 log2(i;^)) , 
oo4 J 

27rr3\/AW^ (^1 + ^(1 + 37;2 + 21og7)+ (8.11) 

+^(15 + 42-^2 + + 4(19 + 9^;2) log 7 + 36(log 7)^ + 12 log2(^;^)) . 
o84 J 

Here we used the asymptotic expansion of the background metric component and the relation 
between and T. We remind that this expansion is valid as long as the velocity satisfies 
e/ilog7 < 1. 
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Figure 9. The ratio of the diffusion constants k±,k\\ to their value in the unflavored background 
are shown as a function of the velocity v, for different values of the perturbative parameter et — 
0,0.25,0.5,0.75 (from bottom up, or from blue to red). This ratio is as well the ratio of the jet 
quenching parameters with respect to their super symmetric values, since they are proportional to the 
diffusion constants. 

The corrections coming from the flavor backreaction in are the same as for the drag 
force with the addition of corrections coming from the worldsheet temperature in equation 
(B.9). We conclude that the jet quenching parameter related to the transverse momentum 
broadening in the flavored background is increased with respect to the value it takes in 
the unflavored background. The same conclusion applies to the momentum acquired in the 
longitudinal direction, since also for gy the corrections are always positive. Visually, this is 
shown in the plots of figure 9, where the ratio of the diffusion constants (equivalently the 
ratio of the jet quenching parameters) is represented. 
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Figure 10. Here, the jet quenching parameters (in units of vJS^T^ ) are shown as a function 

of the velocity v, at a fixed temperature, and for different values of the perturhative parameter th ~ 
0,0.25,0.5,0.75. 

In figure 10 the jet quenching parameters are plotted, relative to the transverse directions 
and the longitudinal direction respectively. 

The Einstein relation in terms of the jet quenching parameter and the drag force can be 
written as 




where we have to use the worldsheet temperature. This relation is well known to be valid in 
the = 4 supergravity dual, as well as in dilaton-gravity models such as [32], and remains 
unchanged in the backreacted background, order by order in eh- 

9 Conclusions 

The physics of quarks propagating through dense media is a subject of intense research 
nowadays. The fact that this problem, most probably, involves QCD at strong coupling 
implies a challenge to the traditional perturbative methods. Holography provides a valuable 
tool, in spite of the fact that the gravity dual model is only qualitatively related to quantum 
chromodynamics. The traditional objections to the presence of supersymmetry are bypassed 
at finite temperature, where supersymmetry is indeed broken. In the present work we have 
analyzed the properties of probes on a background that was constructed some time ago [1]. 
This background is analytic. However, the price to pay is that it is a perturbative deformation 
in powers of a parameter which is roughly the Veneziano parameter ~ Nf/Nc- The use of this 
dual geometry only makes sense in the spirit of effective field theory. The background comes 
with a cutoff r^,, beyond which the predictions blow up because of the presence of a Landau 
pole. Conversely all the physical predictions at scales much lower than the cutoff r <^ r^, 
are robust and unique up to corrections which die off with powers of r/r*. In the present 
context, the lowest cutoff energy that one can consider is given by the plasma temperature 
(or the horizon radius r > r^). When examining the physics in the furthest possible IR, it 



makes sense to discard the UV corrections, or equivalently send — )■ oo, and work with the 
effective metric written in (2.6)-(2.10). To make sense of the analytic expansion in terms of the 
parameter et, related to the Veneziano factor, the effective energy bound is r < ta ~ r/^e^/'^'^ 
This implies that the mass of the quark probes cannot be taken to be very large, namely the 
relevant range for masses is Mq < 3.5 GeV. 

Given all this, and moreover, that among all possible "conifold like" models we have chosen 
the simplest one, namely AdS^ x S5, we shall only trust qualitatively our results. At a 
practical level this amounts to finding the sign of the first order correction to physically 
motivated quantities. Whether it is positive, negative or vanishing is all we can say within 
the present effective approach. Of course, as emphasized in the text, there is not a unique 
prescription as to what two theories should be compared. The choice both here and in [1], [2] 
has been to keep the total amount of degrees of freedom constant, as measured by the energy 
density. Using the entropy density would have given qualitatively equivalent results. The 
phenomenological outcome of this work can be summarized in few sentences. The effective 
thermal masses (constituent Mc and kinetic Mkm) increase in a flavored theory as compared 
with their values in an unflavored one. The meson melting transition occurs at a lower value 
of the ratio Mq/T. For fixed quark mass this implies a higher temperature. The qq screening 
length is also diminished. The conductivity is reduced, which is consistent with the fact 
that the drag force increases. The fiavor diffusion constant is lowered, but the broadening of 
jets, as parametrized by the jet quenching parameter, is enhanced. All in all, the consistent 
picture that arises seems to point out in the direction of saying that the fundamental degrees 
of freedom, have a larger cross section than the adjoint ones. 
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A Derivation of the conductivity 

We will mimic [16] to obtain the conductivity for the probe D7-brane. To simplify notation 
let us write for the pullback metric 

dsl =5«(r)dt2+5,,(r)dxW + (7,,(r,x(r))dr2 + dsL(x(r),n , (A.l) 
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where O"" = (9, ^, ip) represent a 3-sphere wrapped by the probe brane, and the pullback metric 
components are related to the background lOd components given by 

G tp'"^ 

9tt = Gtt , Qxx = Gxx ) drr = Grr H 2~'l ^ ' (^■^) 

dsl = %(1 - i^') [{W'f + + [W'f] + ^IL_^(1 - ^2)2(^3^2 (A.3) 

We have field strength F = dA where, as explained in the main text, the gauge field is given 
hy A = Ao{r)dt + {E^t + ^^(r))dx. 

The action for the probe brane has the following form in Einstein frame 

S = Sdbi + Scs = -T7 j d\ ^-det(58 + e-*/2F) + Scs , (A.4) 

where dsg = ds\ is the induced line element in the Einstein frame, and F = 2iTa'F. When 
evaluating the determinant in the DBI action, the compact space will factorize since F does 
not take values along these directions, the contribution reads (with a,b = 9, "0, <p) 

detgat = GlJl - i^^f (1 + - 1) (1 - ^')) ^ - 5.n.(r, • (A.5) 

In this expression we have separated the radial dependence, gintifjipir)), from the angular 
part. Altogether 

grr{e'^\9tt\glx - glxEl) + \gtt\glx^'x - alx^'S) 9int—^ , (A. 6) 

hence 



Sdbi = -V^M j dtdr ^ {grr\gtt\gxx - e"* {gwEl - \gtt\A';^ + gxxA'^)) glx9int , (A.7) 

where Af = 27r'^N'j:T7, V3 = J d^x and we have used the factorization of the angular variables 
to obtain the volume of the 3-sphere, J d^9 s\n9 /S = 2it'^. Notice the length dimensions 
[A^] = -8 and [gint] = 6. 

There are first integrals of motion for Ao^x , which are holographically related to the charge 
and current densities 

S^DBI S^DBI 



(Jx) , (Ai 



with dimensions [ng] = [jx] = —3, since [Cdbi] = —5. Dividing by 27ra'M and defining the 
quantities hg = ng/(27ra'M) and (Jx) = (Jx) /{27ra'M), we can express the radial derivatives 
of Aq^x as 



A' 
^0 



\gtt\grr e _ \g^^\g^^ 

(A.9) 



'\ gxx \ \gtt\glx9int + e ^ds-^l^q - gxx{Jx)) 



= -f^A', . (A.IO) 
\9tt\ 
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One then defines a special locus (which we name singular shell) as the radius r^j such 
that the numerator in the square root of Aq vanishes, and fixes the value of (Jx) to make 
the denominator vanish at the same point, to keep the gauge field real (and, therefore, the 
Legendre transformed action w.r.t Aq^x stays real too) 



e E^- \gtt\gxx 



\gtt\gxx9int + e {\gtt\nq- gxx{J'. 



, 
. 



(A.ll) 
(A.12) 



Equation (A.ll) defines implicitely the sigular shell r^{rh). Notice that for a metric in the 
form given in (A. 2) we have, from equation (A.ll), that the singular shell and the horizon 
radius are related by 

= + ■ (A.13) 

From equations (A.ll -A.12), and using the definitions given in this appendix, we obtain the 
expression 

27ra'K, . 



{2TTa'M)^gxxgint + e-'^nU ^ 



Which gives the following expression in terms of the bulk components: 



f^DC 



Ex 



(A.14) 



(A.15) 



'(27ra')W2G,.,G3 (1-^2)3 1 + 



^ - 1 ) ) + (27ra')2e-*n2Gr2 



Taking the limit i^^: — is equivalent to taking the limit rg- — s- r^, where 

^2V^^ (^-i)(^-* 



,,2\3 



{27rarM'GxxG'(l-r) 



2\3 
h) ) 



and also 



(27ra')2e-*n2G 



XX 



4e" 



n„ 



^(0)^2^4 



(A.16) 



(A.17) 



which gives the result quoted in (5.2). As a final remark concerning dimensions, notice that 
for the conductivity, [ax] = —1, is consistent with Ohm's law, since [j^] = —3 whereas 
[E] = [F, 
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Following again the procedure stablished by Karch and O'Bannon in [16] we can study the 
drag force from the conductivity. For very heavy carriers we expect a classical equation of 
forced motion 

f = -W + / (A.18) 
In the stationary situation, with f = Ex, this will give rise to a steady current 

jx = (TxEx = UgV, (A. 19) 
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with Uq the density of charge carriers and v their velocity. In the very massive probe hmit 
we can discard the first term in the expression (A. 14) and compare it with(A.19), obtaining 



(A.20) 



Using this result and the definition of the singular shell (A. 11), as well as the explicit values 
of the induced metric components in terms of the bulk metric components we find that the 
velocity and drag force are given by 



\Gtt\ 



f 



1 



27ra' 



(A.21) 



B Drag Force and Jet Quenching 

Central to the discussion of all these effects is the action for a fundamental string 

where the string worldsheet coordinates in the static gauge = t, = r are labelled 
by a,/3 = 0, 1 and Gmn is the bulk metric in the Einstein frame defined in equation (2.2), 
with components given by equation (2.9). From here the momentum currents can be readily 



(B.l) 



derived 



vr?. 



27ra' 



(B.2) 

-9 \GuGxxx' 

Let us consider first the stationary solution: the ansatz for the classical trailing string reads 

X'^=vt + x{r), X' = 0, i = 2,3, (B.3) 

where v is the velocity of the quark moving along the direction and x{r) determines its 
(steady state) position. The classical on-shell action for the trailing string only depends on 
the radial derivative of x. Therefore, the momentum conjugate to x, vr^, is constant along 
the radial direction. Moreover, vr^ is related to the momentum loss of the quark, due to the 
friction of the medium as the quark is moving in the direction, vr^ = —dpi/dt. This is the 
drag force that we want to compute. 

To proceed, we solve the equation for x' in terms of the constant tTx and of the background 
metric. This yields 



x' = 27ra'7rT 



Grr 



Gtt + Gxxv^ 



GttGxx (27ra'7ri)2 + e'^GuGxx 



(B.4) 



Reality of this expression is ensured if the denominator and the numerator both change sign 
simultaneously at the same point defined by 



I Gtt I 



G, 



dpi 
dt 



27ra' 



■e^'^GxxV 



(B.5) 



- 28 - 



From the first equation and (2.6) we can solve for in terms of rh as = ^/^rh ■ with 
7 = (1 — u/2)~^/^. Also, from the second, we have the following implicit expression for the 
drag force 

(B.6) 



^ - dt - 27ra'^ 

Comparing (B.5) and (B.6) with (A. 21), we find consistency provided r„ = rg. 

In order to proceed with the fluctuations it is most convenient to go to a different gauge. 

~^a={a, = t + arUi = r) ; CV) = ^ = ^7^^%^, (B.7) 
for which the induced metric is diagonal for the same trailing string ansatz as (B.3) 

9ab=\ ^ ^ ^ ,2^ [xx'G^^f . (B.8) 

\ Gtt + X^Gxx I 

One sees that is the location of an event horizon of the induced worldsheet metric 5'oo('^s) = 
0, for which, performing the standard manipulations one can compute the associated Hawking 
temperature 



rj. _ J_ _ / ^'tt + ^xx - Git {G'ix + Gtt^) 



/3s 47r V GrrGuGl, 



(B.9) 



from which the result given in (7.3) follows upon expansion up to second order in This, as 
usual, defines the worldsheet temperature Tg in terms of the bulk temperature T. 

In order to compute the jet quenching parameter, we have to resort to the analysis of 
fluctuations. The ansatz is 

X'^ = vt + x{r) + 6X^{t,r) ; = 6X' i = 2,3. (B.IO) 

Expanding (B.l) to second order we have (in (Tq coordinates) 

= / rfV (^g;'6xljxl + st^K^xi^ ' (B-ii) 

where the (diagonal) tensors G^^^^ are given in the main text (8.5) (8.6). Now we can derive 
the equations of motion in this gauge. Let di = d^i 

dr (e^/^gVdr^i) - a;2e*/2g00$, = I =±, \\ (B.12) 

with 

/ X2'3(^o,ai = r) = e--°^x(a;,r) 

I Xi(a°,ai =r) = e^-^°^||(w,r) ' ^ ' 
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or, explicitly, 



H 

were H and Z were given in the main text (8.5) and (8.6). 
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